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We  study  the  thermal  coupling  potential  between  a  high  temperature  metal  hydride  (MH)  tank  and  a 
Solid  Oxide  Fuel  Cell  (SOFC)  aiming  towards  the  design  of  an  efficient  integrated  system,  where  the 
thermal  power  produced  during  normal  SOFC  operation  is  redirected  towards  the  MH  tank  in  order  to 
maintain  H2  desorption  without  the  use  of  external  heating  sources.  Based  on  principles  of  thermody¬ 
namics,  we  calculate  the  energy  balance  in  the  SOFC/MH  system  and  derive  analytical  expressions  for 
both  the  thermal  power  produced  during  SOFC  operation  and  the  corresponding  thermal  power  required 
for  H2  desorption,  as  a  function  of  the  operating  temperature,  efficiency  and  fuel  utilization  ratio  in  the 
SOFC,  and  the  MH  enthalpy  of  desorption  in  the  tank.  Based  on  these  calculations,  we  propose  an  in¬ 
tegrated  SOFC/MH  design  where  heat  is  transferred  primarily  by  radiation  to  the  tank  in  order  to 
maintain  steady-state  desorption  conditions.  We  develop  a  mathematical  model  for  this  particular 
design  that  accounts  for  heat/mass  transfer  and  desorption  kinetics  in  the  tank,  and  solve  for  the  dy¬ 
namics  of  the  system  assuming  MgH2  as  a  storage  material.  Our  results  focus  primarily  on  tank  operating 
conditions,  such  as  pressure,  temperature  and  H2  saturation  profiles  vs  operation  time. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Considerable  research  has  been  invested  in  recent  years  towards 
the  development  of  non-petroleum  energy  carriers  for  use  both  in 
mobile  (transportation,  portable  devices  etc)  and  stationary  appli¬ 
cations  (backup  and  auxiliary  power  units,  distributed  electricity 
generation  etc).  In  this  direction,  hydrogen  has  emerged  as  an 
attractive  option  due  to  its  production  from  a  variety  of  primary 
energy  sources,  its  high  energy  content  and  clean  exhaust  product. 
Furthermore,  the  use  of  H2  can  significantly  contribute  to  the 
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successful  penetration  of  renewable  energy  (RE)  sources  in  the 
electricity  grid,  as  it  offers  a  promising  alternative  medium  for  RE 
storage,  which  could  be  subsequently  used  as  a  clean  fuel  for  ve¬ 
hicles,  as  well  as  for  distributed  electricity  production  through  the 
use  of  fuel  cells,  or  even  internal  combustion  engines,  etc. 

Efficient  H2  storage  remains,  however,  a  significant  technolog¬ 
ical  barrier  towards  the  widespread  application  of  hydrogen  pow¬ 
ered  devices  and  vehicles,  as  the  volumetric  energy  density  of 
uncompressed  hydrogen  gas  is  very  low.  In  this  direction,  metal 
hydrides  (MHs)  are  being  considered  as  promising  hydrogen  stores 
due  to  their  inherent  high  gravimetric  hydrogen  content,  and  sig¬ 
nificant  effort  is  devoted  to  the  development  of  MH-based  upscaled 
storage  systems  and  their  testing  (in  terms  of  storage  capacity, 
charging/discharging  kinetics,  reversibility  and  cycling,  etc.)  at 


AG.  Yiotis  et  al.  /  Journal  of  Power  Sources  26 9  (2014)  440-450 


441 


different  operating  conditions  [1,2  .  As  a  matter  of  fact,  commercial 
MH  tanks  became  very  recently  available  for  stationary  applications 
[3]  building  on  the  distinct  advantages  of  MgH2,  i.e.  its  high 
reversible  gravimetric  capacity  (7.6%  w/w),  flat  reaction  plateau  at 
low  pressures,  and  low  cost  due  to  the  abundance  and  facile 
extraction/processing  of  Mg. 

Hydrogen  storage  in  MH  tanks  is  an  exothermic  process  that 
produces  significant  excess  heat  (in  the  range  of 40—80  kj  mol-1  for 
most  materials  tested  in  the  literature).  Unless  the  produced  excess 
heat  is  removed,  it  results  in  increasing  tank  temperatures  that 
eventually  satisfy  the  activation  energy  barrier  for  H2  desorption, 
and  thus  the  material  can  no  longer  absorb  more  hydrogen  leading 
to  incomplete  tank  chargings.  The  design  of  efficient  cooling  sys¬ 
tems  for  MH  tanks  is  thus  a  crucial  factor  for  exploiting  the 
maximum  theoretical  storage  capacity  of  the  material. 

Hydrogen  desorption,  on  the  other  hand,  is  an  endothermic 
reaction  that  requires  sufficient  heat  fluxes  to  be  provided  to  the 
MH  bed  in  order  to  maintain  the  desired  desorption  mass  fluxes. 
Attempts  to  develop  integrated,  self-sustainable  systems,  where 
the  excess  thermal  power  produced  from  fuel  cells  under  normal 
operating  conditions  is  redirected  towards  an  MH  tank  to  sustain 
H2  desorption  have  been  recently  reported.  Most  of  these  works 
have  focused  primarily  on  the  thermal  coupling  between  low 
temperature  MH  tanks  and  Polymer  Electrolyte  Membrane  (PEM) 
fuel  cells,  where  the  cooling  air  flow  stream  from  the  fuel  cell  is 
used  to  heat  the  MH  tank  [4,5]. 

For  high  temperature  MH  materials,  such  as  MgH2  that  exhibits 
high  thermodynamic  stability  (requiring  -74  kj  moH1  of  desorbed 
H2)  and  high  activation  energy  (desorption  kinetics  are  sufficiently 
fast  at  temperatures  T  >  300  °C),  thermal  coupling  is  only  feasible 
with  high  temperature  fuel  cells,  such  as  Solid  Oxide  Fuel  Cells 
(SOFC)  that  typically  operate  at  temperatures  above  600  °C.  For 
such  a  combination  of  high  temperature  MHs  coupled  with  SOFCs, 
different  heating  solutions  have  been  considered  so  far,  including 
the  after-burning  of  unused  H2,  and  the  electrical  heating  of  the  MH 
tank  [6],  However,  the  design  of  a  self-sustainable,  integrated,  high 
temperature  SOFC/MH  tank  system  that  can  efficiently  maintain  H2 
desorption,  until  the  tank  becomes  practically  fully  depleted,  still 
remains  an  open  scientific  and  engineering  problem. 

In  this  work,  we  propose  such  an  integrated  system  design 
where  a  MgH2  multiple  tank  system  is  heated  through  radiation 
and  conduction  by  a  SOFC  operating  at  a  typical  500  W  electrical 
power.  We  study  the  energy  balance  between  the  SOFC  and  the 
storage  system,  and  derive  analytical  expressions  for  the  thermal 
power  produced  during  SOFC  operation  and  the  corresponding 
thermal  power  required  for  H2  desorption  from  the  tanks,  as  a 
function  of  the  fuel  utilization  rate,  the  SOFC  operating  temperature 
and  efficiency,  and  the  MH  enthalpy  of  desorption.  Based  on  these 
calculations,  we  develop  a  mathematical  model  that  accounts  for 
heat  transfer  in  the  integrated  system,  as  well  as  mass  transfer  and 
desorption  kinetics  in  the  tanks,  and  solve  for  the  dynamics  of  the 
coupled  system  at  typical  operating  conditions  and  different  fuel 
utilization  rates.  Our  results  focus  primarily  on  tank  operating 
conditions,  such  as  tank  pressure,  temperature  and  H2  saturation 
profiles  vs  operation  time. 


room  temperatures,  SOFCs  typically  operate  at  temperatures  above 
600  °C,  where  YSZ  exhibits  increased  ionic  conductivity.  In  fact, 
Joule  heating  at  the  ceramic  electrolyte  is  the  primary  source  of 
energy  dissipation  in  SOFCs,  as  shown  in  Appendix  A. 

Hydrogen  oxidation  occurs  at  the  three-phase  boundary  (TPB: 
porous  space-ionic  conductive  material-electronic  conductive  ma¬ 
terial)  of  the  anode  GDE; 

H2(g)  +  02--H20(g)  +  2e-  (1) 

The  electrons  produced  at  the  TPB  are  then  transferred  over  the 
external  circuit  towards  the  cathode  GDE,  where  oxygen  is  reduced 
to  ions; 

02(g)  +  4e~  — >202~  (2) 

The  ions  are  then  in  turn  transferred  towards  the  anode  through 
the  electrolyte  to  support  the  first  reaction. 

The  overall  reaction  is  strongly  exothermic  with  an  enthalpy 
change  approximately  equal  to  AH/  =  -250  kj  moU1  (standard 
enthalpy  of  water  formation).  However,  given  that  the  entropy  of 
the  system  is  reduced  during  water  formation,  only  part  of  this 
energy  is  available  for  producing  electrical  (non-expansion)  work 
according  to  the  change  of  the  Gibbs  Free  Energy  AG  =  AH  -  TAS 
under  isothermal  conditions. 

From  a  thermodynamic  point  of  view,  the  reaction  energy  that 
can  be  used  to  produce  electrical  work  as  a  function  of  the  reacting 
species  concentration  at  the  anode  and  cathode  TPBs  is  described 
by; 


AG  =  AG0  +  RT\n(  Pli2°?n  )  (3) 

\Pn2,aPoJ 

where  P  =  P/P°  is  the  partial  pressure  of  the  species  normalized 
with  the  standard  pressure  Pq  =  1  atm  that  corresponds  to  AGo.  The 
associated  reversible  Fuel  Cell  potential  is  then  given  by  the  well- 
known  Nernst  equation; 


e'-e»-§H 


H20  ,a 


Ph2Jo'cJ 


(4) 


where  AG  =  -zFE,  z  =  2  is  the  number  of  electron  moles  exchanged 
per  H2  mole,  and  F  is  Faraday’s  constant. 

In  principle,  the  maximum  energy  conversion  predicted  by  Eqs. 
(3)  and  (4)  is  never  achieved  due  to  irreversibilities  in  the  fuel  cell 
associated  with  the  diffusion  and  mixing  of  species  in  the  GDE’s, 
ohmic  resistances  in  the  electrodes  and  the  electrolyte,  and  the 
activation  energies  of  the  reaction  kinetics.  Due  to  these  effects, 
part  of  the  Gibbs  Free  Energy  of  Eq.  (3)  is  irrersibly  converted  to 
heat.  These  irreversibilities  (typically  termed  overpotentials, 
namely,  concentration  Vc,  ohmic  V0,  and  activation  Va,  respectively) 
are  manifested  by  reducing  the  Nernst  potential  Er  (which  corre¬ 
sponds  to  the  maximum  theoretical  voltage  of  the  SOFC)  and  are 
increasing  functions  of  the  current  density  in  the  fuel  cell; 


2.  Solid  Oxide  Fuel  Cell  energy  balance 

SOFCs  are  electrochemical  devices  that  efficiently  convert  the 
Gibbs  Free  Energy  of  hydrogen  oxidation  directly  to  electrical  work. 
They  are  constructed  by  an  ion  conducting,  gas-tight  ceramic 
membrane  (typically  yttria-stabilized  zirconia  (YSZ)),  that  acts  as 
an  electrolyte,  sandwitched  between  two  porous  gas  diffusion 
electrodes  (GDE)  with  both  ionic  and  electronic  conductivity  [7,8], 
Due  to  the  high  electrical  resistivity  of  such  ceramic  materials  at 


E(i)  =Er  -  V0(i)  -  Vc(t)  -  Vfl(i)  (5) 

A  sample  calculation  of  the  overpotentials  based  on  a  mathe¬ 
matical  model  developed  by  Zhu  et  al.  [9]  for  a  typical  2D  planar 
SOFC  cell  is  included  in  Appendix  A.  This  analysis  demonstrates 
that  ohmic  overpotentials  in  the  electrolyte  are  the  primary  source 
of  energy  dissipation  in  SOFCs  operating  at  maximum  electrical 
power.  This  is  an  expected  result,  since  both  the  electric  conduc¬ 
tivity  of  GDEs  is  relatively  high,  and  the  activation/concentration 
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overpotentials  are  negligible  at  maximum  electrical  power  (namely 
at  medium  current  density).  Furthermore,  at  such  elevated  tem¬ 
peratures  the  activation  barrier  of  the  electrochemical  reactions  is 
easily  overcome,  leading  also  to  negligible  activation  overpotential 
effects. 

Based  on  this  simplification,  we  can  obtain  analytical  expres¬ 
sions  for  the  energy  balance  in  an  integrated  SOFC/MH  system  as  a 
function  of  operating  temperatures,  fuel  utilization,  enthalpy  of 
desorption  and  the  temperature  of  the  exhaust  gases  exiting  to  the 
environment.  Denoting  by  ri  the  molar  flow  rate  of  H2  from  the  tank 
to  the  SOFC  anode  and  by  a  the  mole  fraction  of  H2  that  actually 
reacts  (utilization  ratio),  the  chemical  power,  Wch,  produced  in  the 
SOFC  is  expressed  as  follows; 

Wch  =  an^Hf\.  (6) 

A  fraction  rj  of  this  power  corresponds  to  the  produced  electrical 
power  Wei  and  the  remaining  corresponds  to  the  thermal  power 
Qth  that  can  be  used  to  heat  the  MH  tanks; 

Wei  =  ?jari|AHf|,  (7) 

and 

Qth  =  (1  -  v)an\/S.Hf\.  (8) 

Part  of  the  produced  thermal  power  is  typically  used  for  pre¬ 
heating  the  reactants  to  the  reaction  temperature  (namely,  the 
SOFC  stack  operating  temperature).  If  air  is  used  as  the  oxidating 
agent  injected  at  the  cathode  at  stoichiometric  analogy  to  the  H2 
molar  flow  rate,  the  required  power  for  heating  the  reactants  is; 

Qheat  =  ^heati  (-0 

where 

Top 

Qheat  =  J  (Cp,H2 +0.5cpo2 +2CpiN2)dT,  (10) 

To 

and  Top  is  the  fuel  cell  operating  temperature,  cPi;  the  heat  capacity 
under  constant  pressure  of  component  i,  and  To  the  initial  tem¬ 
perature  of  the  reactants  at  the  inlet  stream. 

This  can  be  further  generalized  by  considering  that  H2  is  sup¬ 
plied  at  a  higher  temperature  than  the  rest  of  the  reactants,  equal  to 
the  average  operating  temperature  of  the  MH  tank,  Ttank; 

Top  T0p 

Qheat  =  /  Cp,H2dr+  /  (0.5cPio2 +2cp,N2)dT.  (11) 

T,  a„k  T0 

At  the  same  time,  the  heat  power  required  to  sustain  H2 
desorption  from  the  MH  tank  is; 

Q.des  =  ri|AHdes|,  (12) 

where  AHdes  is  the  enthalpy  of  H2  desorption  at  the  tank  operating 
temperature. 

Besides  the  thermal  power  Qu  ,  produced  at  the  SOFC,  part  of  the 
thermal  power  contained  in  the  outlet  stream  of  the  reaction 
products  (exhaust  gases)  can  be  also  used  for  heating  either  the  MH 
tank,  or  for  preheating  the  reactants  inlet.  Assuming  that  the 
products  of  the  reaction  exit  from  the  SOFC  to  the  environment  at  a 


temperature  equal  to  T\,  the  usable  thermal  power  in  the  exhaust 
gases  stream  is; 

Qexh  =  TlQexh:  (13) 

where 

Top 

Qexh  =  J  [“cp,H20  +  0.5(1  -  a)cp  02  +  (1  -  “)Cp,H2  +  2Cp  N2)]dT. 
T, 

(14) 

Furthermore,  heat  can  be  also  produced  using  part  of  the  elec¬ 
trical  power  Wel  of  the  SOFC  when  all  other  heat  sources  cannot 
satisfy  the  heat  flux  for  H2  desorption  from  the  MH  and  for  pre¬ 
heating  the  reactants  in  the  SOFC.  Assuming  that  a  fraction  y  of  the 
produced  power  is  used  from  joule  heating,  then  the  produced 
thermal  power  is; 

Qj  =  yweh  (15) 

and 

Pd  =  (1  -  y)We,  (16) 

is  the  remaining  usable  electrical  power. 

The  operation  of  the  SOFC  is  thus  sustainable  (reaches  a  steady 
state)  when  its  temperature  remains  constant  at  Top  and  the 
following  relation  is  satisfied; 

Qth  +  Qexh  +  Qj  +  Qenv  Qheat  Qdes  =  (12) 

where  Qenv  is  the  heat  power  emitted  to  the  environment  in  the 
form  of  waste  heat.  In  the  above  notation,  the  SOFC  is  heated  by  an 
external  source  when  Qenv  >  0,  and  cooled  by  the  environment 
when  Qenv  <  0. 

Denoting  qenv  =  Qenv/ri  we  then  take; 

(1  —  “F  ^exh  “F  Qenv  ~  Q heat  —  qdes  —  9-  (13) 

From  Eqs.  6—17  and  given  that  Qth  =  Wch  -  VVe(,  we  have; 

a((Y  —  l)7?  “F  1 ) |  AHj I  —  |AHdes|  +  Qexh  —  qheat  ~  ~Qe nv  (19) 

Net  energy  is  released  to  the  environment  when  Qenv  <  0, 
which  leads  to; 

“((7-  1)4  +  1)|aH/|  -  |AHdes|  +  qexh  -  qheat  >  0  (20) 

Note  that  in  the  above  analysis  no  assumption  has  been  made 
regarding  the  efficiency  of  the  SOFC,  namely  which  fraction  of  the 
useful  chemical  power  is  converted  to  electrical.  Such  an  approxi¬ 
mation  can  be  made  by  taking  into  account  that  both  concentration 
and  activation  overpotentials  are  negligible  compared  to  ohmic 
overpotentials  when  the  SOFC  operates  at  its  maximum  power  We/, 
as  shown  in  Appendix  A.  Under  these  conditions  the  electrical 
power  of  the  SOFC  is; 

We,  =  El  =  (Er  -  V0)I  =  Erl  -  l2r ,  (21 ) 

where  r  is  the  ohmic  resistance  of  the  electrolyte,  /  the  current  that 
passes  through  the  circuit,  Er  the  open  circuit  voltage  of  the  SOFC 
(see  Eq.  (5))  and  VVe(  r  =  Erl  =  an  AGp  is  the  maximum  reversible 
electrical  power  that  can  be  produced  by  the  SOFC,  neglecting  ir¬ 
reversibilities.  The  actual  peak  electrical  power  of  the  SOFC  is  then 
obtained  at; 
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^  =  E,-2Ir  =  0  (22) 

Namely,  the  peak  power  is; 

We/,max  =  §.  (23) 

at  /max  =  £r/2r.  The  SOFC  efficiency  can  be  thus  calculated  at  /max  as 
follows; 


^e/.max  WeimSLX 

AGj  | 

Ef  /Ar 

A  Gf  | 

Wch  Welr 

A  Hf\ 

£2/2r 

A  Hf\ 

A  Gr 
=  0.5 — - 

AH/ 


(24) 


Given  thus  a  typical  SOFC  efficiency  y  =  0.5  AGy  /  Affy-  ,  we  can 
determine  the  lowest  utilization  ratio  of  H2  in  the  SOFC  to  ensure 
the  release  of  net  energy  in  an  integrated  SOFC/MH  system  for 
typical  operating  conditions  and  assuming  MgH2  as  a  hydrogen 
storage  material.  By  combining  Eq.  (20)  with  Eqs.10  and  11,  and 
setting  T0  =  300  K,  W  =  600  K,  Top  =  1000  K,  A Hdes  =  75  kj  mol-1, 
A Gf=  -200  kj  mol-1  (calculated  at  Top),  A Hf=  -250  kj  mol  1  and 
the  heat  capacities  of  H2,  02,  N2,  and  H20(g)  referred  by  NIST  [10], 
we  deduce  that  the  lowest  utilization  ratio  for  net  energy  produc¬ 
tion  must  be  a  >  0.73,  when  none  of  the  electrical  power  at  the 
SOFC  is  used  for  heating  the  MH,  namely  when  y  =  0.  In  the 
opposite  case,  where  all  the  chemical  power  of  H2  oxidation  is 
irreversibly  converted  to  heat  in  the  SOFC,  namely  when  y  =  1  (or 
equivalently  when  =  0),  the  minimum  fuel  utilization  ratio  for  net 
energy  production  is  a  >  0.43. 


Fig.  1.  Schematic  of  the  integrated  SOFC/MH  system  considered  in  this  study.  Only  one 
quarter  of  the  total  system,  that  consists  of  16  MH  tanks  in  total,  is  shown  here.  The 
computational  domain  is  highlighted  with  red  color.  Heat  is  transferred  from  the  SOFC, 
that  operates  at  T  =  1000  K,  to  the  inner  surface  of  the  metal  ring  primarily  by  radi¬ 
ation.  The  ring  serves  to  minimize  temperature  gradients  around  the  tanks  and  provide 
uniform  desorption  patterns,  while  the  air  layer  blocks  heat  transfer  by  conduction,  so 
that  the  tanks  operate  at  significantly  lower  temperature  than  the  SOFC.  (For  inter¬ 
pretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  this  article.) 


3.  Integrated  system  design 

The  above  analysis  demonstrates  that  thermal  coupling  in  in¬ 
tegrated  SOFC/MH  tank  systems  is  possible  in  terms  of  the  available 
thermal  power,  even  in  the  case  of  MHs  with  high  enthalpies  of 
formation  (such  as  MgH2),  provided  a  high  hydrogen  utilization 
ratio  of  a  >  0.73.  Additional  analysis,  however,  is  required  in  order 
to  study  the  dynamics  of  integrated  systems  in  terms  of  heat/mass 
transfer  and  H2  desorption  kinetics,  and  ensure  the  operating 
conditions  where  heat  transfer  from  the  SOFC  to  the  MH  tank  oc¬ 
curs  fast  enough  to  maintain  both  a  constant  H2  feed  to  the  SOFC 
and  steady  state  tank  operating  conditions  (i.e.  practically  constant 
tank  pressure  and  temperature). 

To  this  end,  we  develop  a  mathematical  model  that  describes 
heat  transfer  in  the  form  of  radiation  and  conduction  from  the  SOFC 
towards  the  MH  tank,  and  heat  and  mass  transfer  within  the  MH 
tank  itself  during  H2  desorption  at  constant  demand.  The  kinetics  of 
H2  desorption  are  described  by  a  first-order  reaction,  where  the 
reaction  rate  and  activation  energy  are  obtained  by  a  least-square 
fit  on  experimental  desorption  data  from  the  literature  [11].  The 
mathematical  model  is  then  solved  in  2D  (assuming  significantly 
higher  gradients  in  the  radial  than  the  axial  dimension  of  the  tank) 
using  COMSOL  Multiphysics  4.3  finite  element  software  package  for 
the  coaxial  system  design  shown  in  Fig.  1. 

Fig.  1  shows  only  one  quarter  of  the  entire  system,  which  is  built 
around  a  cylindrical  SOFC,  0.1  m  in  diameter  and  0.1  m  in  length. 
The  SOFC  is  surrounded  by  a  thick  stainless  steel  ring  that  contains 
[16]  equally  spaced  MH  cylindrical  tanks  0.04  m  in  diameter  and 
0.1  m  in  length  each.  The  SOFC  and  the  metal  ring  are  seperated  by 
an  air-filled  region  0.07  m  in  thickness  that  acts  as  a  low  conduc¬ 
tivity  barrier  between  the  SOFC  and  the  MH  tanks.  The  dimensions 
and  operating  conditions  are  based  on  a  commercially  available 
Almus  Ubocell  SOFC  model  capable  of  500  W  electrical  power  [12], 


This  particular  design  is  selected  in  order  to  take  advantage  of 
the  radiant  heat  fluxes  from  the  SOFC  to  the  surrounding  MH 
tanks,  while  maintaining  a  significantly  lower  operating  tem¬ 
perature  for  the  tanks  (typical  operating  temperature  for  the 
tanks  300  °C  <  T  <  400  °C,  while  the  typical  operating  temper¬ 
ature  of  the  SOFC  is  well  above  600  °C)  by  introducing  the  low 
conductivity  air  layer  between  the  SOFC  and  the  metal  ring.  This 
allows  for  reduced  operating  pressures  within  the  tanks  and 
minimizes  aging  of  the  MH  material,  while  assuring  sufficiently 
fast  desorption  kinetics.  Furthermore,  the  metal  ring  surrounding 
the  tanks  allows  for  highly  conductive  heat  transfer  around  the 
tanks  that  produces  uniform  temperature  profiles  at  the  tank 
walls  and  thus  uniform  H2  profiles  within  the  tanks. 


4.  Heat  transfer  between  the  SOFC  and  the  tanks 

The  radiant  heat  flux  from  the  walls  of  the  SOFC  stack  to  the 
internal  perimeter  of  the  metal  ring,  S,  of  the  porposed  integrated 
design  is  a  function  of  temperatures  of  the  SOFC  and  the  metal 
surface  of  the  ring  and  the  emissivity  of  the  materials  [13,14]; 

Efc  ^tank  l  ftank  J 

where  A/c  is  the  external  surface  area  of  the  fuel  cell,  Atank  the  in¬ 
ternal  surface  area  of  the  surrounding  ring  which  faces  the  fuel  cell, 
efc,  t'tank  are  the  emissivity  coefficients  of  the  SOFC  and  tank  sur¬ 
faces,  respectively,  and  a  is  the  Stefan— Boltzmann  constant. 

The  continuity  of  heat  fluxes  at  the  internal  perimeter,  S,  of  the 
ring  then  reads; 
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~n  ■  (ksVT\s+ -  kaVT \s_)  =  qfc  (26) 

where  ks,ka  are  the  metal  and  air  thermal  conductivity,  respectively, 
and  J T  is  the  unit  vector  normal  to  the  surface,  pointing  from  the 
metal  ring  towards  the  air  layer. 

Similarly,  the  conservation  of  heat  within  the  air  layer  and  the 
metal  ring  reads; 

PlCpl^  +  V-(k,VT)=0  (27) 

where  i  =  s,a  ( s  for  steel  and  a  for  air)  and  p;  is  the  density  and  CPi  is 
the  specific  heat  capacity  of  the  medium,  respectively. 

In  the  following  simulations,  we  take  e/c  =  etank  =  0.05, 
ka  =  0.08  W/(m  K),  Cpa  =  1075  J/(kg  K)  for  air,  and  ks  =  24  W/(m  K), 
Cps  =  475  J/(kg  K),  and  ps  =  7700  kg  m  3  for  the  metal  ring.  The 
above  are  average  values  at  T  =  700  K,  while  the  air  density  is 
calculated  by  the  ideal  gas  law  at  1  atm  pressure. 

5.  Metal  hydride  tank  modeling 

In  our  simulations,  we  consider  that  the  cylindrical  tanks  are 
filled  with  MgH2  powder  which  is  a  very  attractive  hydrogen 
storage  medium  due  to  its  high  storage  capacity  that  ideally  may 
reach  up  to  7.6%  w/w,  and  thus  results  in  high  energy  densities. 
MgH2  however  exhibits  significant  thermodynamic  stability  and 
high  activation  energy  for  desorption,  thus  resulting  in  high  tank 
operating  temperatures  (300  °C  <  T  <  400  °C).  This  latter  charac¬ 
teristic  makes  it  an  ideal  candidate  for  coupling  with  high  tem¬ 
perature  SOFCs. 

The  MgH2  powder  desorbs  hydrogen  according  to  the  following 
reaction; 

MgH2(s)— >Mg(s)  +  H2(g),  AHdes  (28) 

where  AHdes  >  0  is  the  heat  of  the  endothermic  reaction  which  is 
-74  kj  mol-1  at  standard  state  conditions. 

The  mass  balance  for  H2  in  the  gaseous  phase  (void  space)  of  the 
MH  bed  reads; 


For  heat  transfer,  we  consider  the  bed  as  a  homogeneous  and 
isotropic  effective  continuum,  where  the  thermal  properties  are 
averaged  over  the  entire  solid  and  void  space  and  convective  heat 
transfer  is  neglected; 

PbCpb|+V-(-/CeVT)  =  mAHdes  (32) 

where  pb  is  the  bed  density,  PMgH2  is  the  density  of  the  metal  hy¬ 
dride,  e°  is  the  initial  bed  porosity  at  maximum  absorption,  Cpb  is 
the  bed  specific  heat  capacity,  I<e  is  the  effective  thermal  conduc¬ 
tivity  of  bed,  and  AHdes  is  the  heat  of  H2  desorption.  The  bed  density 
is  a  function  of  the  absorbed  H2  mass  fraction  /3  =  pHs/PHS,sat  as 

Pb  =  Pmsh2  (1  -  e°)  fl  -  ^  (1  -  0) '  J ,  where 

PHS,sat  =  (^-e0)i^^Ps,sat  is  the  maximum  concentration  of 
absorbed  H2.  It  is  then  straightforward  to  show  that  the  bed 
porosity  is  £  =  1  -PblPs.  where  Ps  =  pMg  +  )3(PMgH2  -  pMg)  is 
approximated  by  a  linear  interpolation  between  the  densities  of 
pure  magnesium  and  pure  MgH2. 

The  enthalpy  of  desorption  is  related  to  the  standard  heat  of 
MgH2  formation,  A H0,  as; 

r 

^^des  —  AF/q  +  J  (Vp  Mg  +  Cp ,H2  —  Cp,MgH2)d7’  (23) 

298K 

where  cp  Mg,  cp  h2,  and  cp ,MgH2,  are  the  specific  molar  heat  capacities 
of  Magnesium,  H2  and  MgH2,  respectively.  Assuming  an  average 
desorption  temperature  of  T  =  300  °C,  we  calculate  an  average 
value  for  the  heat  of  desorption  A Hdes  =  77.9  kj  mol  1  which  is  used 
in  this  study. 

It  should  be  noted  that  the  effective  conductivity  of  the  bed 
depends  strongly  on  the  particle  size  distribution  of  the  metal  hy¬ 
dride  powder,  ranging  from  2  to  8  W  nT1  K_1  for  unmilled  and 
mildly  ball-milled  powders  with  particle  sizes  in  the  order  a  few 
micrometers,  to  0.1—1  W  itT1  K_1  for  nanoparticle  powders  [16,17], 


^  +  V-(pg(u»  =  —  m 
where 

<^>  =  ~VPg 
Pg 


(29) 


(30) 


5.1.  Kinetics  of  H2  desorption 

We  assume  first-order  reactions  kinetics  for  H2  desorption  from 
MgH2  [18,15]; 

m  =  PHS,sat ^  =  -Ko  exp  PpHStsat  (34) 


where  pg  is  the  H2  gaseous  phase  mass  per  unit  bed  volume,  which 
is  related  to  the  H2  density  pg  through  the  porosity  e,  as  pg  =  epg,  m 
is  the  rate  of  H2  desorption  per  unit  bed  volume,  k  is  the  bed 
permeability,  Pg  =  pgRTIMg,  pg  is  the  viscosity  of  H2  with  a  tem¬ 
perature  dependence  calculated  by  Sutherland's  formula 
Mg  =  Po7o  +  C/T  +  C(T/T0)3/2,  C  =  72  K,  To  =  293.85  K, 
Ho  =  8.76  x  10~6  Pa  s,  T  is  the  temperature  and  R  is  the  ideal  gas 
constant.  We  note  here  that  properties  with  '  are  averaged  over  the 
bed  volume,  V,b,  (solid  and  void  space)  and  rfi  =  tiMg/Vbed. 

Assuming  negligible  solubility  and  diffusivity  of  gaseous  H2  in 
the  solid  phase  (MgH2  +  Mg  mixture)  [15],  the  mass  balance  for  H2 
in  the  metal  hydride  reads; 


where  pHS  is  the  volume  averaged  concentration  of  chemisorbed  H2 
in  the  bed,  which  is  related  to  the  concentration  per  unit  solid 
volume  as  pHS  =  (1  -  e)pHS. 


where  Ea  is  the  activation  energy  and  K0  is  a  rate  constant.  We 
should  note  however  that  this  is  merely  a  first  order  approximation 
of  the  desorption  process  that  may  take  place  in  several  steps, 
involving  several  reaction  mechanisms  [18], 

In  order  to  calculate  the  activation  energy  and  the  desorption 
rate  constant,  we  apply  a  least  square  fit  of  the  integrated  Eq.  (34)  to 
the  experimental  data  of  Cabo  et  al.  1 11  ]  who  performed  desorption 
measurements  in  a  wide  temperature  range  of  330-390  °C  and 
fitted  the  respective  results  to  the  Avrami— Erofeev  equation.  These 
results  are  consistent  with  other  literature  data,  i.e.  [19,18].  The 
integrated  linear  equation  is  /n(l-a)  =  -Kt,  where  a  =  1-/3  and 
I<  =  K0exp(-Ea/RT),  where  we  fit  only  data  in  the  range  0.2  <  a  <  0.8. 

Fig.  2  shows  the  fitted  curves  and  the  experimental  points 
where  the  slope  of  the  curves  expresses  I<  =  K(T).  The  activation 
energy  is  found  Ea  =  153.8  kj  mor1  using  again  a  least-square  fit  for 
1C(T)  values  on  the  Arrhenious  equation,  in  very  good  agreement 
with  values  reported  in  the  literature  for  unmilled  and/or  mildly 
ball-milled  powders  [19,18,20], 
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Fig.  2.  Kinetics  of  hydrogen  desorption  at  different  temperatures.  Points  represent 
experimental  data  by  Cabo  et  al.  [11]  and  lines  are  the  least-square  fitted  curves  for 
first  order  reaction  kinetics.  The  slopes  of  the  curves  express  the  desorption  rate 
constants  K  as  a  function  of  temperature. 

Eq.  (34)  is  then  modified  to  account  for  the  effect  of  H2  con¬ 
centration  in  the  void  space,  which  is  known  to  limit  desorption  at 
increased  pressures  [21].  Given  that  hydrogen  absorption  is  a  first- 
order  reaction  at  least  for  T  <  340  °C  and  mild  hydrogen  pressures 
[20,21],  it  is  convenient  to  express  the  overall  desorption/absorp¬ 
tion  rate  per  bed  volume  as; 

fn  =  -K0e(  E‘/RT)f>Hs(^^)  (35) 

where  Ea  is  the  activation  energy  for  H2  desorption  and  Peq  is  the 
equilibrium  pressure  of  H2,  calculated  using  Van't  Hoffs  expression 

[22] 

\n(Peq/P0)  =  A H0/RT  -  A S0/R  (36) 

where  Pq  is  the  atmospheric  pressure  and  ASo  is  the  standard  en¬ 
tropy  change  of  the  system. 

6.  Results  and  discussion 

The  above  mathematical  model  is  solved  using  COMSOL  Multi¬ 
physics  4.3  in  a  computational  domain  that  corresponds  to  1/32  (or 


7r/16  in  cylindrical  coordinates)  of  the  entire  integrated  system 
(shown  in  red  color  in  Fig.  1 )  assuming  continuity  boundaries  in  the 
axial  direction.  We  perform  a  series  of  numerical  simulations  under 
typical  operating  conditions  of  the  integrated  SOFC/MH  system  in 
order  to  study  both  the  dynamics  of  the  H2  desorption  process  and 
the  overall  storage  system  operation  for  different  fuel  utilization 
ratios. 

In  a  typical  simulation,  we  assume  that  the  individual  MH  tanks 
and  the  surrounding  space  are  initially  preheated  at  350  °C,  and  we 
impose  a  constant  temperature  Top  =  727  °C  at  the  perimeter  of  the 
SOFC.  The  external  perimeter  of  the  metal  ring  is  thermally  insu¬ 
lated,  ~n  VT  =  0,  and  continuity  boundary  conditions  are  applied  at 
the  wall  of  each  the  MH  cylinder. 

We  then  solve  the  transient  model  until  t  =  3600  s,  while  there 
is  no  H2  demand  from  the  MH  tanks,  as  we  impose  a  no  flux 
boundary  condition  at  the  wall  of  each  MH  cylinder.  For  times 
t  >  3600  s,  we  impose  a  constant  mass  flux,  N,  at  the  wall  of  each 
MH  cylinder  equal  to  N  =  ri/A  =  Wei/(32a7]nrtLAHf),  that  satisfies 
the  desired  electrical  power  from  the  SOFC.  Here  rt  is  the  tank 
radius,  A  is  the  total  surface  of  the  MH  tank  walls,  and  we  impose 
Wei/L  =  5000W/m  (such  as  the  integrated  system  with  a  total 
length  of  L  =  0.1  m  will  produce  500  W  electrical  power).  Simula¬ 
tions  are  terminated  after  2.5  h  of  tank  operation. 

Fig.  3  shows  the  residual  H2  saturation  profiles  in  a  single  MH 
tank  at  four  different  time  steps  of  applied  constant  mass  flux  at  a 
typical  a  =  0.8  fuel  utilization  in  the  SOFC.  These  profiles  show  that 
regions  closer  to  tank  walls  empty  faster  than  regions  closer  to  the 
tank  center,  as  the  local  temperature  in  the  former  region  is  ex¬ 
pected  to  be  higher  thus  leading  to  faster  desorption  kinetics.  This 
is  also  evident  in  Fig.  4  that  shows  the  temperature  profiles  at  the 
same  time  steps,  rescaled  by  the  maximum  and  minimum  tem¬ 
peratures  at  each  time.  The  corresponding  extreme  temperature 
values  are  also  shown  under  each  profile.  These  profiles  also 
demonstrate  the  effect  of  placing  the  MH  tanks  in  a  highly 
conductive  metal  ring,  that  smooths  out  temperature  gradients  at 
the  perimeter  of  the  tanks,  thus  leading  to  more  or  less  homoge¬ 
neous  desorption  profiles. 

Fig.  5(left)  shows  the  residual  H2  saturation  profiles  in  a  tank  as  a 
function  of  time  for  four  different  utilization  ratios.  It  is  apparent 
that  at  lower  H2  utilization  ratios,  the  MH  tanks  empty  faster  as 
they  need  to  supply  higher  mass  fluxes  to  the  SOFC,  in  order  to 
produce  a  fixed  value  of  electrical  power.  At  the  same  time,  the 
corresponding  heat  fluxes  entering  the  tank  walls, 
Qw/L=  (l/L)fA[rt -V(ksT)]dA,  shown  in  Fig.5-(right),  are  signifi¬ 
cantly  higher  at  lower  utilization  ratios.  Interestingly  enough,  heat 


Fig.  3.  Residual  H2  saturation  profiles  in  a  single  MH  tank  at  four  different  time  steps  of  applied  constant  mass  flux  for  a  =  0.8  fuel  utilization  ratio.  Note  that  the  lower  side  of  the 
tank  is  closer  to  the  SOFC  walls. 
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Fig.  4.  Temperature  profiles  in  a  single  MH  tank  at  different  time  steps  of  applied  constant  mass  flux  for  a  =  0.8  fuel  utilization  ratio.  Note  that  the  lower  side  of  the  tank  is  closer  to 
the  SOFC  walls. 


Fig.  5.  (Left)-Residual  H2  in  the  MH  tanks  vs  time,  and  (Right)-Heat  flux  normal  to  tank  walls  for  different  fuel  utilization  ratios. 


fluxes  appear  to  reach  a  steady  state  value  very  soon  after  the 
application  of  a  constant  mass  flux  boundary  condition  at  the 
perimeter  of  the  tank  under  typical  operating  conditions.  It  is  in  fact 
straightforward  to  show  that  the  steady-state  heat  flow  rate  normal 
to  the  tank  walls  is  equal  to  Qdes  =  riAHdes,  and  by  replacing  with 
the  imposed  mass  flux  at  the  tank  walls  and  assuming  a  typical 
SOFC  efficiency  equal  to  V  =  0-5.  we  find 

'®W/m,  This  corresponds  to  94/aW/m  for 


the  half  MH  tank  of  our  computational  domain,  in  agreement  with 
the  results  shown  in  Fig.  5(right). 

Regarding  the  operating  conditions  of  the  MH  tank  itself,  Fig.  6 
shows  the  tank  pressure  (left)  and  temperature  (right)  for  the  same 
values  of  fuel  utilization.  Assuming  that  a  minimum  value  of  H2 
pressure  of  -1  bar  is  required  to  be  injected  to  the  SOFC  anode  for 
the  proper  operation  of  the  fuel  cell,  it  is  apparent  that  the  tanks 
become  quickly  pressure  depleted  at  lower  fuel  utilization  ratios, 
even  though  there  may  still  be  substantial  quantities  of  absorbed 
H2.  In  the  case  of  a  =  0.75  for  example,  the  tank  pressure  drops 


Fig.  6.  Dynamics  of  tank  pressure  (left)  and  temperature  (right)  at  different  fuel  utilization  ratios. 
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Fig.  7.  (Left)  Residual  H2  saturation  along  the  diameter  of  the  MH  tank  at  different  operation  times  for  a  =  0.8.  Left  side  of  the  plot  is  closer  to  the  SOFC.-(Right)  Porosity  of  the  bed  at 
the  same  times. 


below  the  critical  H2  pressure  threshold  in  less  than  2  h  of  opera¬ 
tion,  while  there  still  remains  0.25  hydrogen  saturation  absorbed  in 
the  tank.  Pressure  depletion  in  the  tanks  at  lower  fuel  utilization 
ratios  (and  thus  higher  H2  mass  fluxes  exiting  the  tanks)  is  asso¬ 
ciated  with  the  poor  conductivity  of  the  MH  bed  (roughly  equal  to 
2  W/(m  K)),  compared  to  the  excellent  conductivity  of  the  sur¬ 
rounding  metal  medium.  Resistance  to  heat  transfer  in  the  tanks 
also  increases  significantly  with  time,  as  the  bed  region  closer  to 
tank  walls  becomes  progressivelly  concentration  depleted  and 
desorption  proceeds  through  bed  regions  closer  to  the  center  of  the 
tank  (see  Fig.  3).  Thus  the  temperature  in  the  central  regions  drops 
in  order  to  satisfy  the  imposed  mass  and  heat  fluxes,  leading  to 
decreasing  tank  pressure,  as  described  by  the  Van't  Hoff  Equation.  It 
is  thus  of  significant  importance  to  maintain  relative  small  tank 
radii,  in  order  to  avoid  early  pressure  depletion  due  to  the  large 
temperature  gradients  within  the  MH  tank. 

The  dynamics  of  H2  desorption  from  the  MH  tanks  are  shown  in 
Fig.  7(left)  as  a  function  of  position  along  the  tank  diameter.  It  is 
apparent  from  this  plot  that  the  H2  desorption  occurs  faster  closer 
to  the  tank  walls,  but  in  a  pretty  much  symmetrical  manner.  This  is 
due  to  the  presence  of  the  highly  conductive  metal  ring  that  allows 
for  practically  uniform  temperatures  around  the  tanks,  regardless 
of  the  wall  distance  from  the  SOFC  which  is  the  heat  source  in  this 
configuration.  This  is  also  apparent  in  Fig.  8  that  shows  the  tem¬ 
perature  in  the  tank  and  the  surrounding  metal  ring  along  the 
diameter  of  the  tank  that  meets  the  center  of  the  SOFC  at  the  same 
time  intervals. 

An  interesting  effect  of  hydrogen  desorption  from  MgH2  is  also 
demonstrated  in  Fig.  7-(right)  that  shows  the  evolution  of  the  bed 
porosity  along  the  tank  diameter  at  the  same  times.  Our  model 


Fig.  8.  Temperature  of  the  metal  ring  and  the  tank  along  the  diameter  at  different 
operation  times  for  a  =  0.8. 


predicts  an  increase  of  bed  porosity  with  time  due  to  the  increasing 
concentration  of  Mg  in  the  bed  mixture,  that  has  a  higher  mass 
density  than  pure  MgH2  (see  also  Table  1). 

The  above  results  verify  that  under  typical  operating  conditions 
a  hydrogen  utilization  ratio  of  a  >  0.8,  representing  the  most 
common  value  in  commercially  available  SOFCs  [23  ],  can  ensure  the 
thermal  coupling  of  SOFC  and  hydrogen  storage  systems,  even  in 
the  worst-case  scenario  of  using  a  metal  hydride  with  a  relatively 
high  desorption  enthalpy,  such  as  MgH2.  However  our  approach 
indicates  that  coupling  can  be  also  achieved  with  even  lower  uti¬ 
lization  ratios  if  one  considers  the  use  of  metal  hydrides  with 
appropriate  desorption  enthalpies.  It  should  be  moreover  noted 
that  although  thermal  coupling  is  mostly  appropriate  for  SOFC 
systems,  the  theoretical  methodology  proposed  in  this  work  can  be 
easily  adapted  to  realistically  describe  a  variety  of  integrated  sys¬ 
tems,  such  as  the  less  demanding  combination  of  low  temperature 
hydride  storage  systems  and  PEM  fuel  cells. 

7.  Conclusions 

In  this  work,  we  study  the  thermal  coupling  potential  between 
a  high  temperature  metal  hydride  tank  and  a  Solid  Oxide  Fuel  Cell 
based  on  principles  of  SOFC  and  MH  thermodynamics.  Civen  the 
high  operating  temperature  of  SOFCs,  we  propose  a  coaxial  system 
design,  where  the  heat  power  of  the  SOFC  is  transferred  primarily 
by  radiation  and  conduction  to  an  assemble  of  MH  tanks,  the 
operating  temperature  of  which  is  reduced  in  the  region 
300-400  °C  (allowing  fast  desorption  kinetics  for  the  case  of 
MgH2,  but  also  several  absorption/desorption  cycles).  Aiming  to 


Table  1 

Physicochemical  properties  of  bed  components  [16,17,22,24]. 


MgH2  and  bed  properties 

Pure  Mg  density,  pMg 

1739  kg  m-3 

Pure  MgH2  density,  ,3Mi!h2 

1450  kg  nrr3 

Mixture  density,  ps 

PMg  +  0(PMgH2  ~  PMg) 

MgH2  molecular  weight,  MMgH2 

26.3  kg  kmol-1 

Specific  heat  capacity 

(1012.6  +  0.483(1  -  273.15)  +  530.97/3)  „4- 

of  bed,  CPb 

Effective  thermal 

2  W/(mK) 

conductivity,  Ke 

Standard  enthalpy  of 

-74.4  kj  mol  1 

formation,  AHo 

Standard  entropy  change,  ASo 

-135.1  k|/(K-mol) 

Bed  permeability,  k 

10- 8  m2 

Initial  bed  porosity,  e° 

0.5 

H2  properties 

Molecular  weight,  Mg 

2  kg  kmol  1 

Density,  pg 

PgMg/RT 

Viscosity,  pg 

Sutherland  Formula 
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Fig.  9.  Schematic  of  the  2D  planar,  cathode-supported  Siemens-Westinghouse  type 
SOFC  used  in  this  study.  Typical  dimensions  for  the  Gas  Diffusion  Electrodes  and  the 
electrolyte  are  taken  by  Weber  at  al  [8], 


calculate  the  thermal  power  required  to  maintain  this  particular 
integrated  SOFC/MH  configuration  self-sustainable,  where  the 
heat  power  produced  at  the  fuel  cell  is  redirected  to  the  MH  tanks 
to  trigger/preserve  H2  desorption,  we  develop  a  mathematical 
model  that  solves  for  both  heat  transfer  from  the  SOFC  to  the  MH 
tanks,  and  heat/mass  transfer  within  the  MH  tanks,  while  ac¬ 
counting  for  H2  desorption  kinetics  from  MgH2.  We  perform  a 
series  of  numerical  simulations  under  typical  operating  conditions 
to  study  the  dynamics  of  desorption  and  the  conditions  under 
which  the  tanks  become  pressure  depleted.  Our  simulations  show 
that  thermal  coupling  is  possible,  even  for  the  case  of  MHs  with 
very  high  enthalpies  of  desorption,  provided  efficiently  high  fuel 
utilization  ratios  at  the  SOFC.  For  poorer  utilization  values,  the 
tanks  soon  become  pressure  depleted  due  to  the  low  thermal 
conductivity  of  the  MH  bed  and  the  resulting  sharp  decrease  in 
operating  temperature,  although  substantial  quantities  of  H2  still 
remain  chemisorbed  in  the  MH  bed.  The  proposed  approach  has  a 
universal  nature  as  it  can  be  also  adapted  for  the  coupling  of  e.g. 
low  temperature  hydrides  with  PEM  fuel  cells. 


v'7h2  =  V-(ph2*u  +7h2)  =th2  (37) 

where  J  Hj  is  the  mass  flux  of  H2  due  to  both  convection  and 
diffusion  in  the  electrode,  U  =  J^Pi^i/p  and  p  =  J2Pi is  the  Darcy 
velocity  andjfensity  of  the  gas  mixture  (where  i  =  H2.H20), 
respectively,  j  H2  is  the  diffusive  mass  flux,  and  rHl  is  the  rate  of  H2 
consumption  at  the  TPB  per  unit  volume  of  the  GDE. 

Assuming  that  the  gas  mixture  is  stagnant,  7?  =  0,  and  the 
effective  diffusivity  of  H2,  Deff,  in  the  mixture  is  related  to  the  mo¬ 
lecular  binary  diffusivity,  Dm,  as  Deff  =  e/rDm,  where  e  =  0.4  is  the 
porosity  and  t  =  1.4  is  the  tortuosity  of  the  GDE,  then  mass  con¬ 
servation  Eq.  (37)  takes  the  following  form; 

V-  (pA?ffVwH2)  =  fn2  (38) 

where  wH2  =  Ph2/p  is  the  mass  fraction  of  hydrogen. 

The  rate  of  hydrogen  consumption  at  the  anode  can  be  related  to 
the  local  current  density  ia  within  the  electrode  using  Faraday's 
law; 


_  SqMh2  . 
r»2  “  2Tla 


(39) 


where  Sa  specific  surface  of  the  TPB  at  the  anode,  MH2  is  the 
hydrogen  molecular  weight  and  F  is  the  Faraday  constant. 

Similar  expressions  hold  for  the  mass  conservation  of  species  at 
the  cathode  GDE,  where  humified  air  is  supplied  at  constant  mass 
flux  at  the  external  surface  and  diffuses  through  the  electrode, 
where  02  is  reduced  to  anions  as  described  by  Equation  (2). 
Equation  (37)  is  thus  valid  for  the  mass  conservation  of  O2,  but 
diffusive  fluxes  are  described  by  the  generalized  Stefan— Maxwell 
expressions  for  a  ternary  system  (namely,  H2O,  O2  and  N2)  [24,25], 
Thus  the  diffusive  mass  flux  of  O2  is; 
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Appendix  A.  Modeling  SOFC  overpotential  curves  and 
efficiency 


N 

j  o2  =  Po2  22  Do2,iVx,-  (40) 

1=1 

where  x,-  is  the  molar  fraction  of  each  species  (i  =  O2.N2.H2O)  and 
B02 ;  are  the  effective  Stefan— Maxwell  diffusion  coefficients  (see 
Ref.  [24]  for  ternary  system  expressions)  by  accounting  the  effects 
of  porosity  and  tortuosity  in  the  molecular  diffusivities,  as  above. 
Note  that  thermal,  pressure  and  forced  diffusion  effects  are 
neglected  in  the  above  expression  [26], 

The  Darcy  velocity  of  the  mixture  can  be  obtained  with  respect 
to  the  N2  diffusive  flux,  given  that  J  N2  =  0; 


For  completeness,  we  provide  the  basic  theoretical  model  for 
the  production  of  electrical  power  in  a  simplified  2D  planar  SOFC 
design  in  order  to  support  our  argument  that  ohmic  resistances  in 
the  ceramic  electrolyte  membrane  is  the  primary  source  of  energy 
dissipation  in  SOFCs  operating  at  peak  electrical  power  resulting  in 
efficiency  »;  =  0.5. 

We  consider  a  cathode-supported  Siemens-Westinghouse  type 
fuel  cell  that  consists  of  an  ion-conducting  (but  gas-tight)  ceramic 
layer,  sandwitched  between  two  porous  GDEs  with  both  electronic 
and  ionic  conductivity  (see  Fig.  9).  Pure  hydrogen  is  supplied  by  the 
MH  tank  at  a  constant  mass  flux  towards  the  external  surface  of  the 
anode  GDE  and  it  diffuses  through  the  porous  electrode  towards 
the  TPB  where  it  reacts  with  O2  to  produce  electrons  and  vapor 
water  (see  Eq.  ( 1 )).  The  produced  water  vapors  then  counter-diffuse 
towards  the  external  surface  of  the  anode  GDE.  By  treating  the 
electrode  as  an  effective  continuum,  mass  conservation  of  H2  at 
steady  state  conditions  reads; 


N 

J  n2  =  Pn2  u  +  j  n2  =  0=>"u  =  -  J]  D jv^-Vx,-  (41) 

i=  1 

The  conservation  of  O2  then  reads; 

V-(po2"u +7o2)  =ro2  (42) 

where  the  Darcy  velocity  in  the  cathode  is  calculated  by  Eq.  (41 ), 
the  diffusive  mass  flux  by  Eq.  (40),  and  the  rate  of  oxygen  reduction 
can  be  expressed  in  terms  of  the  cathode  local  current  density,  ic,  as 
follows; 


ro2 


scMq2 

4  F 


C 


(43) 


where  Sc  specific  surface  of  the  TPB  at  the  cathode  and  Mq2  is  the 
oxygen  molecular  weight. 


Table  2 

Variables  for  Butler— Volmer  equations. 
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Anode  electrode 

Pre-exponential  factor,  AHz 

Activation  energy,  JSjjg8 

Site  density,  F 

Sticking  coefficient,  S^2 

Pre-factor  for  exchange  current  density,  kHz 
Activation  energy  for  exchange  current  density,  £h2 
Cathode  electrode 

5.59  x  1015  m2/(s  mol) 
88.12  kj  mol"1 

2.6  x  1CT5  mol  m_1 

0.01 

2.07  x  109  Am"2 

87.8  kj  mor1 

Pre-exponential  factor,  A0z 

Activation  energy,  EjJ*s 

Pre-factor  for  exchange  current  density,  ko2 
Activation  energy  for  exchange  current  density,  Eq2 

4.9  x  1013  Pa 

h2 

5.19  x  108  A  ITT2 

88.6  kj  moG1 

The  current  density  produced  at  the  TPB  of  both  the  anode  and 
cathode  GDEs  is  typically  expressed  as  a  function  of  the  local  acti¬ 
vation  overpotentials  Va  at  the  electrolyte/electrode  interface  using 
the  Buttler— Volmer  equation  [27,28],  In  this  contribution,  we  adopt 
the  following  formalism  proposed  by  Zhu  et  al.  [9,29]; 

For  the  anode  electrode  (hydrogen  oxidation); 


ia  —  io  ,n 


exp 


-  exp 


(-0.5FVa.aJ 


(44) 


where 


!o  ,a 


=  k. 


h2  exp 


(ph2/ph2) 

|1/4(Ph2o/p°)3/4 

1  +  1 

S/PH2) 

1 

and 


(45) 


io,c  =  l<02  exp 


eo2) 

1 

(PO; 

>/Pb2) 

.1/4 

RT  ) 

1  +  1 

(po2/po2) 

(48) 


and 


/  £des\ 

Pb2  =^o2  expf  --gr-J  (49) 

The  definitions  and  values  for  the  variables  of  the  above  equa¬ 
tions  are  provided  in  Table  2. 

Assuming  a  spatially  homogeneous  TPB  in  the  GDEs,  a  charge 
balance  for  electrons  in  the  electron  conductive  side  of  the  interface 
then  reads; 

V’  (  —  Ke,effV<pc)  =  Si  (50) 

where  cpe  is  the  electronic  potential,  /ce,eff  is  the  effective  electric 
conductivity,  S  is  the  specific  surface  area  of  the  GDE  equal  to  Sa  for 
the  anode  and  Sc  for  the  cathode,  and  i  is  the  current  density  due  to 
the  reaction  at  the  TPB  as  calculated  by  Eq.  (44)  and  Eq.  (47)  for  the 
anode  and  the  cathode,  respectively. 

A  similar  charge  balance  for  the  oxygen  ions  in  the  electrolyte 
side  of  the  interface  reads; 

V-(-*q,effV«)  =Si  (51) 

where  q>i  is  the  ionic  potential  and  /qeff  is  the  effective  ionic  con¬ 
ductivity  of  the  GDE. 

A  charge  balance  in  the  electrodes  (outside  the  GDEs)  reads; 
V-(— /ceV<pe)  =  0  (52) 


,  _^h2T  ^ 2ttRTMH2 
P h2  ~ - cO - exP 

n2 


and  for  the  electrolyte; 

(46) 

V-(-KiV<Pi)  =  0 


(53) 


The  definitions  and  values  for  the  variables  of  the  above  equa¬ 
tions  are  provided  in  Table  2. 

For  the  cathode  electrode  (oxygen  reduction); 


ic 


!o,c 


exp 


^0.5fVa/j 


-  exp 


^-0.5FVa,cy 


where 


(47) 


where  ne  and  «,■  is  the  charge  conductivity  in  the  electrodes  and  the 
electrolyte,  respectively. 

Eqs.  (37)— (53)  are  solved  subject  to  appropriate  boundary 
conditions  using  a  numerical  model  developed  for  COMSOL  Mul¬ 
tiphysics.  Here,  we  take  fixed  mass  fractions  u>  for  all  species  at  the 
external  surface  of  the  anode  and  cathode  GDE's,  and  fixed  elec¬ 
trical  potential  E  at  the  external  boundaries  of  the  metal  electrodes, 
in  order  to  solve  for  current  production  within  GDE's  and  the 


Fig.  10.  Effect  of  operating  temperature  and  H2  molar  fraction  at  the  anode  inlet  on  the  electrical  power  (open  symbols)  and  voltage  (closed  symbols)  produced  at  the  SOFC  as  a 
function  of  current  density. 
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species  fluxes  across  the  external  CDE's  boundaries.  We  use  the 
potential  difference  E  across  the  electrodes  as  an  independent 
variable  in  order  to  calculate  the  corresponding  current  density 
through  the  external  circuit  and  the  produced  electrical  power. 

The  overall  current  density  produced  by  the  SOFC  is  then 
calculated  at  the  anode  GDE  by  /  =  ^  fjadv,  where  A  is  the  surface 
of  the  anode  GDE  open  to  the  H2  inlet.  The  produced  electrical 
power  is  calculated  as  VVe(  =  El  and  the  fuell  cell  efficient  is 
77  =  Wel/fi\Hf,  where  ri  is  the  calculated  mass  flux  of  H2  entering 
the  anode  GDE. 

Fig.  10  shows  typical  polarization  curves  and  the  produced 
electrical  power  at  different  SOFC  operating  temperatures.  Inter¬ 
estingly  enough,  the  electrical  power  curves  correspond  to  a  second 
degree  polynomial  of  current  density  in  agreement  with  Eq.  (22), 
where  the  internal  resistance  of  the  SOFC  perfectly  matches  that  of 
the  ceramic  electrolyte  layer  of  the  same  thickness  as  a  function  of 
temperature.  Namely,  the  electrical  conductivity  of  the  electrolyte 
lq  =  (ff/T)exp(-E(/(i?T)),  where  <7  =  3.6xl05  S  K  citT1  and 
Ei  =  8  x  104J  mol-1  [9],  The  calculated  fuel  cell  efficiency  at  peak 
electrical  power  is  77  =  0.5  at  all  temperatures. 
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